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ABSTRACT 

It has been previously proposed that some distant open clusters in the Milky Way may have been accreted 
during a dwarf galaxy merger, perhaps associated with the same event that led to the formation of the Galac- 
tic anticenter stellar structure (GASS), also known as the "Monoceros Ring". We have obtained VI and 
Washington+DD(95 1 photometric and medium resolution (R ~ 8000) multi-fiber spectroscopic data for the 
three distant old open clusters Berkeley 29, Saurer 1, and vdB-Hagen 176 (BH 176). These clusters are spa- 
tially coincident with GASS, but radial velocities and spectroscopic metallicities had not been available during 
previous studies of the GASS candidate cluster system. Similar data for the clusters Berkeley 20 and Berke- 
ley 39 have been obtained for calibration purposes. We provide the first reliable radial velocity for BH 176 
(V he , io = 1 1 .2 ± 5.3 km s" 1 ). We also find that V heHo = +95.4 ±3.6 and +28.4 ± 3.6 km s" 1 , for Saurer 1(A) and 
Berkeley 29, respectively. We show that a-enhanced isochrones, while spectroscopically motivated, provide a 
poor fit to Be29 in contrast to previous findings. We find that the clusters Berkeley 29 and Saurer 1 are con- 
sistent with the previously reported characteristics for GASS candidate clusters and the GASS stellar stream 
as derived from M-giant observations. However, the radial velocity and photometric metallicity ([Fe/H] ~ 0.0 
dex) for BH 176 suggests that a connection of this cluster with the putative GASS cluster system is unlikely. 
We reassess the age-metallicity relation for the most likely members of the GASS clusters system for which 
spectroscopic metallicities are now available. 

Subject headings: Galaxy: structure - Galaxy: disk - galaxies: interactions - Galaxy: open clusters and asso- 
ciations: individual (Berkeley 20, Berkeley 29, Berkeley 39, Saurer 1, BH 176) - Galaxy: 
globular clusters: individual (BH 176) 



1. INTRODUCTION 

Frinchaboy et al. (2004, hereafter F04) noted that the out- 
ermost open clusters in the Milky Way seem to lie in a string- 
like configuration that can be fit to an orbital plane coinci- 
dent with the Galactic anticenter stellar structure (GASS). The 
presence of this arc-like structure was inferred from excesses 
of various types of stars (Newberg et al. 2002; Yanny et al. 
2003; Ibata et al. 2003; Majewski et al. 2003; Crane et al. 
2003; Rocha-Pinto et al. 2003, 2005) beyond the apparent 
limit of the Galactic disk, and has also been used to argue for 
the presence of a distinct, extended stellar structure wrapping 
around the disk at low latitudes. The feature has been referred 
to as the "Monoceros Ring" (Newberg et al. 2002; Yanny et 
al. 2003) or GASS (Majewski et al. 2003; Rocha-Pinto et al. 
2003; Crane et al. 2003). However, because of its unfortu- 
nate location behind considerable extinction, the system's true 
shape, orientation, extent, etc. have been difficult to ascertain. 
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Even the location of the structure's center (presumably corre- 
sponding to a "nucleus") remains uncertain and controversial. 
For example, this stellar stream has been argued to be associ- 
ated with the postulated Canis Major (Ibata et al. 2003; Martin 
et al. 2004) or Argo dwarf galaxies (Rocha-Pinto et al. 2005). 

Previous work (Crane et al. 2003) on this stellar arc has 
determined some characteristics including: (1) a velocity- 
longitude trend that indicates a slightly non-circular orbit, 
(2) a velocity dispersion that is smaller than even that of 
disk stars, and (3) a wide metallicity spread from [Fe/H]= 
-1 .6 ± 0.3 dex (Y03) to at least [Fe/H]= -0.4 ± 0.3 dex. F04 
found that at least five globular clusters and one open cluster 
(Tombaugh 2) have positions and radial velocities (RVs) sug- 
gesting an association with GASS (see also Martin et al. 2004; 
Frinchaboy et al. 2005), but pointed to other clusters with spa- 
tial coincidence with GASS that still have unknown RVs. As 
a result, we sought to derive the kinematics and chemistry for 
these other clusters. In this work, the clusters Berkeley 29 
(Be29), Saurer 1 (Sal) 8 , and vdB-Hagen 176 (BH176) from 
the F04 list of postulated GASS clusters are investigated. Ta- 
ble 1 gives the positional data for these clusters, as well as 
for old, outer open clusters Berkeley 20 (Be20) and Berke- 
ley 39 (Be39), which we also observed as metallicity calibra- 
tors. Medium resolution spectra for all five clusters are used 
to determine their bulk RVs and metallicities via measure- 
ment of the Ca II infrared triplet lines. These data are used 
to test whether these clusters are consistent with the observed 
dynamical and chemical trends of the GASS system. While 
Yong et al. (2005) and Penarrubia et al. (2005) suggest the 

8 Originally described as Saurer A in Frinchaboy & Phelps (2002) and 
Carraro & Baume (2003). 



2 



Frinchaboy et al. 



TABLE 1 
Cluster Basic Parameters 



Cluster 


"2000.0'' 


^OOO.o" 


ICY 


b(°f 


Diam." 


E(V-I) puh 


Ref 


vdB-Hagen 176 (BH176) 


15:39:05.4 


-50:03:01.7 


328.4100 


+ 4.3418 


2.3' 


0.70 


1 


Berkeley 20 (Be20) 


05:32:37.0 


+00:11:30.0 


203.4803 


-17.3711 


2.0' 


0.16 


2 


Berkeley 29 (Be29) 


06:53:04.2 


+16:55:39.0 


197.9493 


+ 7.9802 


2.0' 


0.10 


3 


Berkeley 39 (Be39) 


07:46:51.0 


-04:40:30.0 


223.5465 


+10.0915 


7.0' 


0.11 


4 


Saurer 1 (Sal) 


07:20:56.0 


+01:48:29.0 


214.6894 


+ 7.3862 


1.3' 


0.18 


5 



References. — (1) Phelps & Schick (2003); (2) MacMinn et al. (1994); (3) Tosi et al. (2004); (4) Carraro et al. (1994); (5) CBVMPR 
'Fundamental data from Dias et al. (2002) 



metallicity and RV data are insufficient to prove membership 
in GASS, they are sufficient to identify clusters with proper- 
ties inconsistent with those observed for GASS. 

In §2 we analyze photometry for the clusters Be20, Be29, 
Be39, Sal, and BH176. This photometry is used to select tar- 
gets for spectroscopy. In §3, we discuss the spectroscopic ob- 
servations and use RVs derived from the spectra to determine 
the clusters' bulk RVs. The derived mean cluster RVs are used 
to clarify the positions of the giant branches of the clusters in 
the color-magnitude diagram and compare our work to pre- 
vious studies in §4. Finally in §5, we explore each cluster's 
likelihood of belonging to GASS and discuss the implications 
these clusters have for GASS and open cluster studies in the 
Galaxy in general. 

2. PHOTOMETRY OF DISTANT OPEN CLUSTERS 

The clusters Be29 and Sal were imaged in the Cousins V 
and / filters with the University of Virginia's 40 inch telescope 
at Fan Mountain Observatory (FMO). We also make use of 
VI photometry of BH176, taken with the LCO DuPont 2.5m; 
these data are described more fully in Phelps & Schick (2003). 
The details of the VI observations are listed in Table 2. The 
Fan Mountain SITe CCD has 2048 x 2048 pixels that are 24 
/im square, resulting in a pixel scale of 0.367 arcsec pixel -1 
and a field of view of 12.67 arcmin on a side. Be20, Be 39 and 
Sal were also observed using Washington (M, T2)+DD05l 
filters 9 with the FMO 40 inch and with the 40 inch Swope 
telescope at the Las Campanas Observatory (LCO; see Table 
3). The Swope CCD has 2048 x 2048 pixels that are 24 /jm 
square, yielding a pixel scale of 0.697 arcsec pixel" 1 , for a 
field of view of 23.7 arcmin on a side. All exposures were ob- 
tained under photometric conditions. During each run, stan- 
dard stars (VI from Landolt 1992, Washington+DD(95 1 from 
Geisler 1996) were observed at various airmasses between air- 
masses 1 and 2 for photometric calibration of the data (see 
§2.2 below). 

The imaging data have been reduced with the IRAF 10 CC- 
DRED package, following standard techniques as described 
in the IRAF CCDPROC documentation. Instrumental mag- 
nitudes were measured with the DAOPHOT and ALLSTAR 
packages, using the point-spread function method (Stetson 
1987). These magnitudes were transformed to the Lan- 
dolt (1992) standard-star system (for V and / observations) 
and to the Geisler (1996) standard system (for data in the 

9 The T2 and / filters used for this study were identical; hereafter we will 
refer to these filters as I. 

10 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 
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FIG. 1. — Dereddened Be20 Washington+DD05 1 photometry color- 
magnitude and color-color diagrams used for spectroscopic target selection. 
Stars selected as open cluster giant candidates from the Washington+DD05 1 
color-color diagram are denoted with stars (if). 



Washington+DD(95 1 system). The data for BH176 are taken 
from Phelps & Schick (2003) and the errors are presented 
there. For the clusters Be20 and Be39, observed only with 
Washington+DD(95 1 photometry, transformation equations 
from Majewski et al. (2000) were used to estimate V and 
V -I. Figures 1-3 show the M, M-I color-magnitude dia- 
gram (CMD) and M-DD05 1 , M-I color-magnitude diagram 
(CMD) for the clusters Be20, Be39 and Sal . 

2. 1 . Selection of Candidates for Spectroscopic Follow-up 

Washington M, I, and DD05 1 photometry can be used to 
discriminate between late-type dwarf and giant stars (Majew- 
ski et al. 2000) by calibrated (M-I,M-DD05l) color-color 
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TABLE 2 
VI Photometry Observing Runs 



Cluster 


Telescope 


Pixel Scale 


UT date 


Exposures 






(arcsec pixeL 


') 


V I 


BH176 

Be29 

Sal 


LCO 2.5-m 
FMO 1.0-m 
FMO 1.0-m 


0.259 
0.367 
0.367 


1998 Jun 30 

2003 Nov 21 

2004 Nov 21 


a a 

1 x 900* 1 x 900* 
I x 900s 1 x 600s 


a Photometry from Phelps & 


Schick (2003). 








TABLE 3 

WASHINGTON+DD051 PHOTOMETRY OBSERVING RUNS 


Cluster 


Telescope 


Pixel Scale 


UT date 


Exposures 






(arcsec pixel -1 ) 


M 


/ DD051 


Be2() 
Be39 
Sal 


LCO 1.0-m 
LCO 1.0-m 
FMO 1.0-m 


0.697 
0.697 
0.367 


2000 Dec 05 1 X 120.V 
2000 Dec 05 1 X 120.V 
2004 Feb 28 1 X 900s 


1 X 120s 2 X 420s 
1 x 120.« 2 x 420s 
1 X 600.« 2 X 1800s 
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FIG. 2. — Same as Figure 5 for Sal. 



FIG. 3. — Same as Figure 5 for Be39. 



diagrams. However to apply this technique, an estimate of the 
foreground extinction is necessary to correct the photometry 
to standard magnitudes. Literature values (presented in Table 
1) have been adopted to implement this color-color based sep- 
aration technique. We transform E(V—I) into E(B-V) using 
the relation derived by Dean et al. (1978). Finally, Majew- 
ski et al. (2000) provides relations to obtain E(M-DD05l), 



E(M-I) and A(M) based on E(B-V) and the ratio R v = 3.1 
for the selective to total absorption. 

Again, Figures 1-3 show the derived (M,M-I) CMDs with 
their respective (M-I,M-DD05l) o color-color diagrams for 
the clusters Be20, Be39, and Sal. The Washington color- 
color "giant region" is selected to separate giant stars from 
dwarf stars based on their surface gravity, as described in Ma- 
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TABLE 4 

Summary of RV Standard Observations 



Star/Sp type 


Fiber 


TDR 


V r ( km s ) 


(T\/ 
V 


HD 18884 


014 


39.82 


-26.50 


1.28 


(M1.5IIIa) 


025 


39.04 


-26.27 


1.47 




028 


30.71 


-27.14 


1.19 




045 


35.45 


-26.38 


1.38 




074 


37.98 


-26.96 


1.44 




086 


26.87 


-26.29 


1.51 




099 


41.07 


-26.33 


1.44 




103 


34.91 


-26.23 


1.42 




118 


33.42 


-26.15 


1.20 






Average 


-26.47 


0.34 






IAU Value 


-25.8 


0.1 






Difference 


-0.7 




HD 150798 


008 


46.36 


-3.41 


1.05 


(K2II-IH) 


019 


47.99 


-3.04 


0.96 




033 


35.22 


-2.77 


1.25 




074 


40.52 


-5.10 


1.63 




091 


40.48 


-3.21 


0.91 




099 


42.19 


-3.10 


1.13 




127 


34.52 


-3.08 


1.23 






Average 


-3.39 


0.78 






IAU Value 


-3.7 


0.2 






Difference 


+0.3 




HD 157457 


020 


43.81 


17.78 


0.81 


(G8III) 


034 


34.79 


17.05 


0.66 




044 


35.51 


16.82 


0.73 




045 


34.56 


17.50 


1.37 




116 


34.56 


17.88 


0.61 




125 


30.61 


18.43 


1.17 




135 a 


11.41 


22.46 


2.09 






Average 


17.48 


0.32 






IAU Value 


17.4 


0.2 






Difference 


+0.1 




HD 161096 


030 


42.29 


-12.29 


1.04 


(K2III) 


038 


35.49 


-12.29 


0.98 




039 


37.83 


-12.45 


1.07 




040 


38.85 


-11.56 


0.97 




048 


31.59 


-12.45 


1.00 




057 


28.02 


-12.41 


0.75 




077 


41.44 


-12.21 


1.12 




098 


43.22 


-12.94 


0.88 




104 


42.25 


-12.69 


0.76 




110 


39.67 


-12.03 


1.05 






Average 


-12.33 


0.37 






IAU value 


-12.0 


0.1 






Difference 


-0.3 





"Due to poorer than average S/N this spectrum was not included in these calculations. 



jewski et al. (2000). For this project, we define the following 
boundaries of the "giant region" for dereddened photometry: 

(M-D) > -0.3 (M-/) + 0.35 for 1.0 < (Af-7) < 1.5 
(M-D)>-0.15 forl.5 < (M-7) < 2.0 (1) 

Stars that fall within the "giant box", are within 4' of the clus- 
ter center (7' for the larger Be39), and are brighter than the 
cluster main sequence turn off (MSTO) are denoted by stars 
(*)• 

3. RADIAL VELOCITIES OF DISTANT OPEN CLUSTERS 

3.1. Data Collection and Reduction 

Spectroscopic data for cluster stars were collected on UT 
2004 March 4-5 using the Hydra multi-fiber spectrograph on 
the 4-meter Blanco telescope at Cerro Tololo Intra- American 
Observatory. The spectra cover A = 7800-8800A with an in- 
strumental spectral resolution of R ~ 8000, or 1.2 A per res- 
olution element. To provide RV calibration, four RV stan- 
dards were also observed, where each "observation" of an RV 



standard entails sending the light down 7-10 different fibers, 
yielding many dozens of individual spectra of RV standards 
(see Table 4). For wavelength calibration, the "Penray" (HeN- 
eArXe) comparison lamp was observed for every fiber setup. 

To get spectra of as many members as possible in these 
compact star systems (diam. ~ 2', see Table 1), we se- 
lected as primary targets cluster red giant branch (RGB) can- 
didates from the cluster CMD. For clusters with Washing- 
ton M, I, and DD051 photometry, these RGB candidates 
were selected based on the CMD and their location in the 
(M-I,M-DD05l) o color-color diagram, as described in §2.1 
(i.e., the stars marked with star symbols in Figures 1- 
3). These stars are the primary targets for spectroscopic ob- 
servation, with stars within 2 ' of the cluster center (R c i ust er) 
having the highest priority for observation. Since the Hy- 
dra spectrograph has "fiber packing" limitations (i.e., there 
is a minimum finer spacing of 2.5"), some fibers were placed 
on "extra" red giant branch candidates (ERGB) selected from 
2' < Rciuster < 4' (see Table 5), and to fill yet more fibers we 
also observed many non-RGB-selected stars in the surround- 
ing field in a search for additional cluster members. Also, 
red clump (RC) and horizontal branch (HB) candidates were 
observed in Be20 and BH176. Spectroscopic targets for Be 
29 and BH176 (clusters for which we were unable to obtain 
Washington photometry) were selected using the CMD of the 
inner 4' of the cluster, with the RC stars for BH176 having 
the highest priority. We also "filled" fibers for these clusters 
in the manner described above. The small numbers of stars 
observed in some fields results from a combination of poor 
observing conditions, which yielded low S/N in the observa- 
tions for some stars, as well as the limitation of the number of 
Hydra fibers that can be packed in a small area. 

Preliminary processing of the two-dimensional images of 
the fiber spectra was undertaken using standard IRAF tech- 
niques as described in the IRAF CCDRED documentation. 
After completing the basic processing, the spectra were re- 
duced using standard routines as listed in the DOHYDRA 
documentation. 

3.2. Radial Velocity Determination 

All radial velocities (V r ) were derived using IRAF's FX- 
COR package, which was first used to determine RVs for the 
standard stars. The RVs of standard stars were measured by 
cross-correlating each Fourier-filtered standard star spectrum 
against every other. The resulting RVs for each individual 
spectra were averaged and the standard deviation measured; 
the results are presented in Table 4. As may be seen, the de- 
rived RVs for the standard stars are all within 1 km s" 1 of the 
IAU values. 

We determine our RV errors from the prescription described 
in Vogt et al. (1995), which is based on analysis of repeated 
standard star spectra (Table 4). The Tonry-Davis Ratio (Tonry 
& Davis 1979, TDR) for each spectrum, measured from FX- 
COR, scales approximately with S/N, which allows us to de- 
termine the radial velocity error using the following equation: 



Error (V r ) ■■ 



a 



(1 + TDR) 



(2) 



where a is a constant calibrated by the standard star data. 
Average velocity uncertainties for the individual Hydra stan- 
dard spectrum are a v ~ 1 km s" 1 for standard star spectra 
with TDR > 25, which is S/N > 20. The cumulative x 2 
statistic for the repeat standard spectrum RV measures is 
S(l +TDR/) 2 (V,,,— {V r }) 2 /a 2 , where (V,) refers to the mean 
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TABLE 5 

Spectroscopy Observing Statistics 



Cluster 


Exposures 


# Observed" 


# RV Stars b 


UTdate 




(total/RGB) 


(S/N > 3) 




BH176 


5 X 1800.S 


34/10° 


19/ T 


2004 Mar. 4 


Be20 


4x 1800.V 


38/12 


20/ 4 


2004 Mar. 5 


Be29 


4 X 1800s 


41/13 


25/ 9 


2004 Mar. 4 


Be39 


5 x 1800.5 


68/39 


54/36 


2004 Mar. 5 


Sal 


4x 1800s 


24/12 


11/ 4 


2004 Mar. 4 



"Total number of stars observed spectroscopically divided into two samples. 1) All stars. 2) Those stars selected as candidate 
RGB members from the cluster CMD (see §2.3). 
b Number of stars with sufficient S/N to obtain an RV measurement. 
c For BH176, "RGB" includes both RGB and RC stars. 




Be20— (g) 



Be20" 




-50 50 100 

V r (km s-') 



50 100 

V r (km s" 1 ) 



FIG. 4. — RV histograms for cluster stars observed, panels (a-e) show 
photometrically selected RGB and RC stars from Washington+DD05 1 color- 
color selection (b,d,e) or from CMD selection in the cluster core (a,c), with 
shaded area denote selected members. Panels (f-j) show RVs for all stars 
observed used to find "new" member stars. NOTE: histogram bins are 5 km 
s" 1 wide. 
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velocity and V r j the velocity of the standard star correspond- 
ing to the z'th observation. We determine a by finding the 
value of x 2 where the probability of exceeding it by chance 
is 50% (defined as Xso)< which approximates la errors. The 
standard star spectra (Table 4) provide a data set with 29 de- 
grees of freedom. From standard tables, we find that xlo = 
28.336 for 29 degrees of freedom, which yields a = 25.4. 

Analysis of the cluster star spectra were also measured with 
FXCOR, and errors were determined from the measured TDR 
using the a derived above. The photometric data and helio- 
centric RV (V r ) and error (Error V r ) measurements for cluster 
target stars with a TDR greater than 5 (S/N > 3) are presented 
in Tables 6-10, along with information on the target selection 
described below in §3.3. The measured cluster star RVs are 
analyzed to determine membership in each cluster. 



FIG. 5. — Dereddened Be20 Washington+Z)ZX)5 1 photometry color- 
magnitude and color-color diagrams. Stars selected as giant candidates from 
the Washington+DD05 1 color-color diagram are denoted with stars Gi- 
ant candidates also observed spectroscopically are denoted with large open 
circles (o) if they are not RV members and filled circles (•) if they are mem- 
bers. Other "fill" stars observed are denoted with open triangles (A) if they 
are not RV members and filled triangles (A) if RV members. 



4. RESULTS 

The photometrically identified, most likely RGB and RC 
stars unambiguously reveal the bulk RV of each cluster, as 
shown by the RV distributions in the left panels (a-e) of Fig- 
ure 4. Since these RGB and RC stars are the most likely to 
provide the "purest" samples of cluster members, we use these 
samples to define one measure of the mean cluster RV and 
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TABLE 6 

Stellar Radial Velocities — BH176 field 





Star 


"2000.0 


c 

c>2000.0 


V 


(V-I) 


TDR 


V T 


Error v r 




Type" 1 




Cluster 














(km s ') 


(km s ') 




Member? 


140 


1 5-38-46 54 


-50-02-00 9 


18.95 


1.840 


10.5 


8.4 


2.2 


fill 


Y 




1 V3R-4fi Q9 


— sn-rn- 1 1 n 

— JU.l/J. J l.W 


19 81 


1 734 


6 5 


—36 5 


3 4 


fill 


N 


1793 


1 V*R-S9 ^4 

1 J.JO.JZ.Jt 


JU.iAIi 1 J .U 


18 08 


2 126 


111 


—1 5 


2 1 


fill 


N 


3313 


ij.jo.j / .tz 


— SOfYVOI Q 
ju.uj.ui .y 


18 99 


1 847 


6.5 


28.4 


3.4 


RC 


N 


3332 


1 J .JO. J /.Jl 




18 90 


977 


17 1 


-40.0 


1.4 


HB 


N 


3994 


1 V^R-SQ 47 


JVi.UJ. JV7.0 


18 80 


1.862 


6.5 


12.6 


3.4 


RC 


Y 


4190 


1 J .J7.UU.UJ 


JU.Vrti 1U. J 


17 83 


2.234 


24.4 


6.4 


1.0 


RGB 


Y 


4245 


15:39:00.24 


-50:01:33.1 


19.16 


0.904 


24.4 


15.7 


1.0 


HB 


Y 


4360 


15:39:00.57 


-50:02:17.0 


18.94 


1.850 


8.4 


9.9 


2.7 


RC 


Y 


5278 


15:39:02.96 


-50:00:14.5 


18.91 


1.831 


5.7 


30.4 


3.8 


ERGB 


N 


5988 


15:39:04.69 


-50:03:17.1 


19.49 


0.998 


5.7 


2.2 


3.8 


HB 


Y 


6194 


15:39:05.14 


-50:05:05.6 


19.25 


0.943 


14.9 


-77.5 


1.6 


HB 


N 


8251 


15:39:10.10 


-50:02:24.0 


18.84 


1.793 


9.2 


-0.5 


2.5 


RC 


N 


8468 


15:39:10.47 


-50:06:33.5 


18.96 


1.874 


6.1 


-6.8 


3.6 


ERGB 


N 


9312 


15:39:12.70 


-50:01:59.9 


18.05 


2.243 


20.2 


11.4 


1.2 


RGB 


Y 


9327 


15:39:12.72 


-50:03:49.6 


18.99 


1.841 


6.9 


16.9 


3.2 


RC 


Y 


10823 


15:39:16.77 


-50:02:20.4 


18.93 


1.901 


7.8 


19.0 


2.9 


ERGB 


Y 


14629 


15:39:29.82 


-50:00:41.2 


18.33 


1.906 


20.2 


-80.8 


1.2 


fill 


N 


15302 


15:39:32.18 


-49:59:51.1 


17.53 


2.319 


10.5 


34.4 


2 2 


fill 


N 



Spectroscopic candidates as selected from the CMD. 
"RGB" are stars along the RGB within the cluster radius. 
"ERGB" are the same only selected from outside the cluster radius. 
"HB" stars are Horizontal branch candidates. 
"RC" stars are red clump candidates. 

"Fill" stars are randomly selected stars from out side the cluster radius. 



TABLE 7 

Stellar Radial Velocities — Berkeley 20 field 



Star 


«2000.0 


<52000.0 


M 


(M-I) 


(M-DD051) 


TDR 


V r 

(km s _1 ) 


Error V r 
(kms- 1 ) 


Type 


Member? 


10308 


5:32:13.32 


0:10:19.1 


17.53 


1.213 


+0.063 


6.9 


1.5 


3.2 


fill 


N 


10617 


5:32:33.19 


0:09:35.0 


17.32 


1.466 


-0.061 


9.2 


0.1 


2.5 


ERGB 


N 


10645 


5:32:34.72 


0:15:38.9 


17.33 


1.385 


-0.048 


10.0 


24.4 


2.3 


fill 


N 


10741 


5:32:36.99 


0:12:00.8 


17.69 


1.456 


-0.015 


6.5 


74.5 


3.4 


RGB 


Y 


10745 


5:32:37.66 


0:17:04.8 


17.95 


1.259 


-0.007 


5.7 


36.5 


3.8 


fill 


N 


10770 


5:32:37.93 


0:11:09.8 


15.12 


1.909 


-0.064 


35.3 


75.5 


0.7 


RGB 


ya 


10794 


5:32:39.16 


0:17:28.7 


17.74 


1.324 


-0.009 


7.8 


77.9 


2.9 


fill 


Y 


10806 


5:32:38.98 


0:09:27.5 


17.08 


0.943 


+0.032 


8.1 


35.3 


2.8 


HB 


N 


10810 


5:32:39.20 


0:10:31.3 


16.00 


1.503 


+0.006 


20.2 


73.7 


1.2 


RGB 


Y 


10850 


5:32:41.03 


0:06:19.4 


18.49 


0.947 


+0.026 


5.7 


17.7 


3.8 


fill 


N 


10851 


5:32:41.53 


0:10:02.9 


16.25 


1.522 


-0.002 


15.9 


79.3 


1.5 


RGB 


Y 


10864 


5:32:42.31 


0:12:29.7 


17.01 


1.050 


+0.032 


5.7 


30.8 


3.8 


HB 


N 


10940 


5:32:46.55 


0:08:42.2 


17.90 


1.405 


-0.006 


6.3 


40.2 


3.5 


ERGB 


N 


10979 


5:32:49.99 


0:11:54.4 


16.80 


1.154 


+0.013 


13.9 


32.5 


1.7 


HB 


N 


10980 


5:32:50.43 


0:16:12.7 


15.46 


1.736 


-0.027 


30.8 


53.8 


0.8 


fill 


N 


11054 


5:32:54.76 


0:17:30.6 


15.48 


1.723 


-0.254 


30.8 


52.6 


0.8 


fill 


N 


11058 


5:32:54.62 


0:10:49.5 


17.62 


1.519 


-0.053 


10.0 


56.6 


2.3 


ERGB 


N 


11094 


5:32:57.52 


0:15:22.7 


17.09 


1.428 


-0.019 


7.2 


31.2 


3.1 


fill 


N 


11158 


5:33:01.22 


0:09:07.3 


18.19 


1.295 


-0.023 


6.1 


-24.3 


3.6 


fill 


N 


11285 


5:33:10.66 


0:07:53.7 


17.92 


1.430 


-0.057 


7.5 


56.5 


3.0 


fill 


N 



"Star in common with Yong et al. (2005). 



observed dispersion (the top values of Table 1 1 with subscript 
"RGB"). 

With this version of the mean RV established for each clus- 
ter, we can search for other potential cluster members iden- 
tified on the basis of similar RV. The RV histograms for all 
stars in the cluster fields with spectra are shown in the right 
panels (f-i) of Figure 4. When the left and right panels of 
Figure 4 are compared we see the usefulness of photometric 
pre-screening to identify the best cluster RGB candidates for 
these sparse, highly field star contaminated clusters: e.g., in 
the cases of BH176, Be20, and Sal the true RV of the cluster 



would remain somewhat uncertain were only the RV distribu- 
tions in the right panels available. We can now rederive bulk 
cluster RVs using all stars with RVs indicating likely cluster 
membership. For this purpose, we limit the cluster member 
range to ± 7 km s" 1 about the mean defined by V r .RCB, val- 
ues guided by the typical individual stellar RV errors and the 
small expected velocity dispersions of these systems (~ 1-3 
km s" 1 ). Because of this imposed selection bias, the disper- 
sions given in Table 1 1 should be interpreted only as a guide 
to the observed tightness of the RVs used to define the mean 
RV, not as the true intrinsic dispersions of the clusters. The 
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TABLE 8 

Stellar Radial Velocities — Berkeley 29 field 



Star 


C»2000.0 


52000.0 


V 


(V-I) 


TDR 


Vr 

(kms- 1 ) 


Error V r 
(km s-> ) 


Type 


Member? 


7m 

zUl 


O.JZ.J J.4U 


io.4y. iz. j 


1 ^ 71 

ID. / J 


1 nfii 


1 < Q 


70 7 

jy.z 


1 c 
l.J 


mi 


IN 


70S 
zuo 


^■s7-07 07 


1 7 7 

io.4y. i / . / 


1 7 in 

1 /JU 


1 1 ah 

1 . 14U 


Q 9 

y.z 


1 Q 1 


7 S 
Z.J 


fill 
1111 


IN 


246 


6-52-50 76 

U.JZ 'W. / V) 


16:49:39.2 


17.28 


0.991 


7.8 


121.2 


2.9 


fill 


N 


249 


6-53-14 68 


16:49:40.7 


17.47 


1.075 


8.1 


-0.5 


2.8 


fill 


N 


462 


6:53:14.92 


16:52:01.3 


17.23 


0.952 


11.1 


61.6 


2.1 


fill 


N 


482 


6:52:50.68 


16:52:17.5 


16.41 


0.960 


13.9 


-27.7 


1.7 


fill 


N 


523 


6:53:25.94 


16:52:38.9 


17.67 


1.021 


6.1 


49.2 


3.6 


fill 


N 


549 


6:53:33.24 


16:52:57.4 


16.94 


0.946 


10.0 


28.5 


2.3 


fill 


Y? a 


725 


6:53:32.80 


16:54:26.1 


16.43 


1.236 


22.1 


31.4 


1.1 


fill 


Y 


822 


fvS^'fll 47 


1U.JJ.U1.7 


16 30 


1 121 


20 2 


27 4 


1.2 


RGB 


Y b 


ooo 


D.jj.Uj.O / 


ID. J J. 1 J.y 


1 A f*(* 


1 . Joo 


i^ 1 


77 ^ 
Z / . J 


n 7 
u. / 


l\ W 1 ) 


I 


Q7fl 
/U 


O.J 5. 15 Ay 


lO.JJ.lO.O 


17 7 1 
1 /.zl 


1 1 71 
l.lzj 


111 
1 1. 1 


£7 7 
/.z 


7 1 

z.l 


fill 
till 


IN 


yj4 


O.J J.UZ.ZU 


ID. J J. JJ. J 




I .U / / 


OA 7 
ZU.Z 


117 
— 11. z 


1 7 
l.z 


KKJD 


IN 


948 


6:53:08.06 


16:55:40.7 


16.68 


1.058 


9.2 


29.1 


2.5 


RGB 




996 


6:53:06.50 


16:55:53.1 


17.49 


1.066 


7.5 


32.5 


3.0 


RGB 


Y 


1000 


6:53:04.37 


16:55:54.3 


14.37 


1.781 


30.8 


26.4 


0.8 


RGB 




1078 


6:52:59.74 


16:56:17.5 


18.00 


1.014 


7.5 


20.8 


3.0 


RGB 


Y 


1165 


6:53:01.67 


16:56:45.3 


15.97 


1.078 


17.1 


53.9 


1.4 


RGB 


N 


1294 


6:53:31.24 


16:57:51.2 


16.57 


1.056 


17.1 


8.4 


1.4 


fill 


N 


1295 


6:53:21.37 


16:57:52.3 


15.81 


0.990 


17.1 


47.6 


1.4 


fill 


N 


1324 


6:53:26.45 


16:58:07.1 


16.30 


1.095 


15.9 


4.2 


1.5 


fill 


N 


1426 


6:52:49.66 


16:58:57.0 


15.93 


1.182 


30.8 


25.8 


0.8 


fill 


Y 


1433 


6:53:25.24 


16:59:00.6 


17.93 


0.860 


6.7 


-17.7 


3.3 


fill 


N 


1437 


6:53:12.25 


16:59:01.3 


16.77 


1.126 


18.5 


24.7 


1.3 


fill 


Y 


1461 


6:53:05.25 


16:59:13.2 


15.31 


1.093 


35.3 


33.2 


0.7 


fill 


Y 



a Star with RV consistent with membership, thought due to CMD loction probably a field star 
^Star in common with Bragaglia et al. (2005). 
c Star in common with Yong et al. (2005). 
d Star in common with CBVMPR. 



mean RVs for the clusters using these enlarged samples are 
also given in Table 1 1 (middle values with subscript "all"). 
Note that while these RV values average over a larger sample 
of stars and therefore might be expected to give a more accu- 
rate measure of the mean cluster RV, it is also the case that a 
higher fraction of interloping field stars with similar RVs may 
enter, and perturb, the average. 

The RVs we have calculated for the five clusters using the 
"RGB" and the "all" samples agree to within 2 km s" 1 , which 
is less than half the observed dispersions. These RVs also 
agree in general with other recent determinations for these 
clusters (bottom values, denoted as "pub" in Table 1 1), but an 
offset of unknown origin is also found for Sal (see §4.4 be- 
low). The CMDs and color-color diagrams for these clusters 
with observed and member stars marked are shown in Figures 
5-8. 

4.1. Isochrone Matching 

We match isochrones to our clusters while constraining the 
metallicity to spectroscopically determined values. We ex- 
plore a range in metallicities centered on the spectroscopi- 
cally determined value for each cluster. Additionally, our ra- 
dial velocity data help clarify the actual positions of the the 
cluster RGBs to improve model isochrone fits for each clus- 
ter. Through these techniques, we re-derive age and distance 
estimates for the clusters Be29, BH 176, and Sal. 

The derived CMDs for these three GASS candidate clus- 
ters are compared to theoretical isochrones from Girardi et al. 
(2000). We use as free parameters the cluster's age, E(V-I), 
(m-M) v , and [Fe/H]. We must stress that given the confused 
appearance of the CMDs for these sparse clusters in heavily 
contaminated backgrounds, especially for Sal, we must be 



careful not to over-interpret the results of these matches. Nev- 
ertheless, useful information can still be extracted from these 
comparisons, and it is a useful exercise to compare these re- 
sults with previously published values. In the analysis that 
follows, visual matching of isochrones covering a range in 
ages and metallicities is carried out after simultaneously vary- 
ing E(V-I) and the distance modulus (m-M) v to get the best 
fit for any given combination of age and metallicity. We have 
used the following equation from Carraro et al. (1999): 

[Fe/H] = log 10 (Z)+ 1.72125 (3) 

to convert between the observed metallicities and the to- 
tal metallicity Z values used for the Girardi et al. (2000) 
isochrones. Additionally, we tested Salasnich et al. (2000) 
a-enhanced isochrone for Be29, see §4.3 below. 

4.2. BH176 

We measure the first reliable RV for BH176 (V r , a u = 1 1.2 ± 
5.3 km s" 1 ) based on 9 stars. This measurement is signifi- 
cantly different and better than our single star measurement 
given in F04, which was of poorer quality and lower resolu- 
tion than our new data. 

The photometric data used for BH176 were obtained from 
Phelps & Schick (2003), where these authors carried out 
isochrone matching using Bertelli et al. (1994) isochrones. 
We rederive the cluster parameters for the same photomet- 
ric dataset, but now use the isochrones from Girardi et al. 
(2000). Figure 9 shows the matched isochrones when we 
adopt Z = 0.008 ([Fe/H] = -0.38), Z = 0.019 ([Fe/H] = +0.0) 
and Z = 0.03 ([Fe/H] = +0.20), with clusters members, deter- 
mined from the spectroscopy, marked as solid triangles. The 
best results are achieved for the highest of these metallicities, 
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TABLE 9 

Stellar Radial Velocities — Berkeley 39 field 



Star 


"2000.0 


52000.0 


M 


(M-l) 


(M-DD051) 


TDR 


V r 

(kms- 1 ) 


Error V r 
(kms- 1 ) 


Type 


Member? 


20289 


7:46:12.10 


^1:39:05.3 


15.26 


1.549 


-0.024 


27.2 


54.7 


0.9 


fill 


Y 


20332 


7:46:13.37 


-4:41:52.3 


17.35 


0.988 


+0.024 


7.2 


34.7 


3.1 


fill 


N 


20362 


7:46:14.97 


-4:37:32.2 


17.32 


1.071 


+0.044 


5.9 


68.0 


3.7 


fill 


N 


20429 


7:46:16.55 


-4:44:59.1 


17.29 


1.098 


+0.037 


8.1 


79.1 


2.8 


fill 


N 


20552 


7:46:21.67 


-4:35:06.1 


17.13 


1.001 


+0.034 


5.5 


55.2 


3.9 


fill 


Y 


20685 


7:46:25.30 


^1:40:33.9 


17.17 


1.305 


-0.032 


11.1 


58.8 


2.1 


fill 


Y 


20729 


7:46:27.25 


^4:32:06.5 


17.28 


1.056 


+0.031 


6.1 


34.4 


3.6 


RGB 


N 


20797 


7:46:28.59 


^1:38:35.4 


17.56 


0.947 


+0.047 


5.5 


37.1 


3.9 


fill 


N 


20814 


7:46:28.29 


-4:46:47.8 


17.73 


1.081 


+0.042 


5.5 


43.8 


3.9 


RGB 


N 


20899 


7:46:30.39 


-4:48:49.6 


16.49 


1 .408 


+0.019 


11.1 


-0.9 


2.1 


fill 


N 


20904 


7:46:32.17 


-4:33:17.0 


17.21 


1.121 


-0.005 


6.1 


-3.9 


3.6 


fill 


N 


21054 


7:46:34.82 


-4:41:14.3 


16.81 


1.407 


-0.039 


14.9 


53.5 


1.6 


RGB 


Y 


21087 


7:46:35.61 


-4:42:09.7 


17.31 


1.081 


+0.028 


6.1 


61.3 


3.6 


RGB 


N 


21091 


7:46:35.88 


-4:40:38.9 


16.90 


1.057 


+0.016 


11.7 


53.4 


2.0 


RGB 


Y 


21152 


7:46:37.18 


-4:40:11.4 


14.33 


1.626 


-0.022 


35.3 


54.8 


0.7 


RGB 


Y 


21206 


7:46:38.58 


-4:39:21.6 


16.49 


1.421 


+0.007 


18.5 


54.0 


1.3 


RGB 


Y 


21247 


7:46:38.99 


-4:41:56.1 


15.09 


1.561 


-0.037 


30.8 


51.5 


0.8 


RGB 


Y 


21305 


7:46:40.53 


-4:38:17.0 


16.62 


1.373 


+0.029 


13.9 


97.7 


1.7 


RGB 


N 


21322 


7:46:39.98 


-4:45:51.3 


17.29 


1.235 


+0.034 


7.5 


67.7 


3.0 


fill 


N 


21365 


7:46:41.27 


-4:40:56.9 


14.65 


1.620 


-0.023 


22.1 


54.8 


1.1 


RGB 


Y 


21366 


7:46:41.45 


-4:39:08.3 


14.66 


1.688 


-0.017 


35.3 


50.9 


0.7 


RGB 


Y 


21426 


7:46:42.49 


^1:41:28.1 


16.95 


1.054 


+0.030 


7.2 


59.4 


3.1 


RGB 


Y 


21432 


7:46:43.30 


^1:34:28.7 


15.89 


1.531 


-0.138 


30.8 


16.0 


0.8 


fill 


N 


21471 


7:46:43.90 


-4:38:13.9 


16.98 


0.938 


+0.035 


9.2 


43.6 


2.5 


RGB 


N 


21523 


7:46:44.57 


-4:40:21.9 


16.73 


1.405 


-0.014 


22.1 


45.8 


1.1 


RGB 


N 


21530 


7:46:44.52 


-4:41:33.3 


17.00 


1.208 


+0.024 


5.3 


55.5 


4.0 


RGB 


Y 


21574 


7:46:45.05 


-4:42:55.0 


17.05 


1.436 


-0.007 


12.4 


50.1 


1.9 


RGB 


Y 


21589 


7:46:45.18 


-4:44:30.2 


17.62 


1.091 


+0.015 


8.1 


45.3 


2.8 


RGB 


N 


21606 


7:46:46.04 


-4:38:58.4 


16.57 


1.423 


-0.029 


27.2 


54.4 


0.9 


RGB 


Y 


21689 


7:46:47.35 


^1:37:08.8 


17.14 


0.997 


-0.016 


5.2 


42.6 


4.1 


RGB 


N 


21695 


7:46:47.99 


^1:32:35.3 


17.72 


1.041 


+0.033 


6.1 


-15.4 


3.6 


fill 


N 


21769 


7:46:48.56 


^4:39:09.1 


16.96 


1.401 


-0.019 


10.5 


50.5 


2.2 


RGB 


Y 


21792 


7:46:48.43 


-4:44:46.8 


16.32 


1.450 


+0.002 


15.9 


37.6 


1.5 


RGB 


N 


21815 


7:46:49.13 


-4:41:55.5 


15.65 


1.540 


-0.029 


24.4 


55.8 


1.0 


RGB 


Y 


21912 


7:46:50.77 


-4:41:28.7 


14.43 


1.584 


-0.017 


35.3 


58.9 


0.7 


RGB 


Y 


21925 


7:46:51.15 


-4:40:29.1 


15.03 


1.721 


-0.038 


35.3 


58.2 


0.7 


RGB 


Y 


21977 


7:46:51.78 


^1:43:18.1 


16.55 


1.522 


-0.004 


22.1 


86.4 


1.1 


RGB 


N 


21982 


7:46:51.74 


^1:44:16.5 


16.40 


1.599 


-0.180 


27.2 


3.7 


0.9 


RGB 


N 


21985 


7:46:52.21 


-4:40:03.9 


17.18 


1.254 


+0.007 


5.9 


36.0 


3.7 


RGB 


N 


22013 


7:46:52.51 


^1:41:14.5 


14.71 


1.588 


-0.023 


35.3 


53.9 


0.7 


RGB 


Y 


22135 


7:46:54.00 


-4:47:51.6 


17.36 


1.170 


+0.063 


7.8 


27.4 


2.9 


fill 


N 


22138 


7:46:53.93 


-4:48:52.8 


16.47 


1.368 


+0.004 


9.6 


-35.3 


2.4 


fill 


N 


OO 1 A A 


/ Ao.jj.Zi 


— 4.JO.J 1 . 1 


1 A 1 Q 

lo. lo 


1 QQO 


— U.UU3 




JO. 3 


1. / 




"V 

I 


22154 


7:46:55.11 


-4:39:27.2 


15.44 


1.563 


-0.026 


30.8 


57.3 


0.8 


RGB 


Y 


22252 


7:46:58.01 


-4:31:40.8 


17.14 


1.045 


+0.023 


6.7 


170.5 


3.3 


fill 


N 


22519 


7:47:02.87 


-4:40:58.8 


17.09 


1.132 


+0.031 


7.8 


72.7 


2.9 


RGB 


N 


22584 


7:47:04.94 


^1:35:42.3 


15.83 


1.509 


+0.002 


27.2 


53.8 


0.9 


fill 


Y 


22603 


7:47:04.88 


^1:40:35.2 


17.40 


1.061 


+0.027 


10.0 


75.4 


2.3 


RGB 


N 


22704 


7:47:06.82 


-4:48:46.7 


17.15 


1.183 


+0.064 


8.1 


56.5 


2.8 


fill 


Y 


22869 


7:47:12.28 


-4:43:17.8 


17.27 


1.195 


+0.014 


7.2 


61.9 


3.1 


fill 


N 


22943 


7:47:15.64 


^1:37:28.7 


17.48 


1.190 


+0.028 


8.8 


75.8 


2.6 


fill 


N 


22989 


7:47:16.75 


^4:39:15.0 


14.51 


1.558 


-0.004 


30.8 


55.8 


0.8 


fill 


Y 


23040 


7:47:17.81 


-4:44:14.4 


14.87 


1.823 


-0.072 


35.3 


62.5 


0.7 


fill 


N 


23082 


7:47:19.41 


-4:42:17.9 


17.43 


1.088 


+0.003 


5.3 


67.1 


4.0 


fill 


N 



which is also the most metal rich value available from Gi- 
rardi et al. (2000). A range of ages (5.6 and 7.1 Gyrs) seems 
to match fairly well the MSTO, the slope of the RGB, and 
the position of the RC, suggesting the age is not well con- 
strained photometrically. Outside this age range a poor match 
is obtained when trying to match these three CMD features 
simultaneously. An E(V—I) from 0.64 to 0.68 is therefore 
derived, resulting in a heliocentric distance of 15.2-15.8 kpc 
(Table 12). We find similar cluster parameters (based on our 
new isochrone matches with updated stellar evolution mod- 
els) to those obtained by Phelps & Schick (2003) using the 
same photometry. 



We agree with the findings of Forbes et al. (2004) that, 
based on current, photometrically-determined [Fe/H], BH176 
does not seem to follow the GASS AMR, in fact it is quite the 
opposite. We have also noted that the RV for BH176 is also 
not strongly correlated to the GASS RV trend, so, collectively, 
the present evidence makes BH 176 unlikely to be a member 
of the putative GASS cluster system. Analysis of the RV find 
that it consistent with an object (open cluster) rotating with 
the flat rotation curve of the Galactic disk. However, this re- 
sult could also be consistent with a globular cluster passing 
through the disk. If BH176 is found spectroscopically to have 
[Fe/H] <~ 0.0, and it is actually a globular cluster (as it has 
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TABLE 10 

Stellar Radial Velocities — Saurer 1 field 



Star 


a 2000.0 


02000.0 


M 


(M-I) 


/Hit nn/ici \ 

(M — DDUDl) 


TDR 


V r 

(kms ') 


Error V r 
(km s ') 


Type 




Member? 


122 


7:20:50.84 


1:43:08.1 


16.89 


1.456 


-0.089 


11.7 


11.7 


2.0 


fill 


N 


145 


7:21:13.56 


1:43:20.6 


16.64 


1.500 


-0.050 


5.3 


88.7 


4.0 


fill 


N 


289 


7:21:11.57 


1:44:49.9 


16.33 


1.484 


-0.114 


12.4 


-6.0 


1.9 


fill 


N 


290 


7:20:46.03 


1:44:52.3 


17.07 


1.601 


-0.195 


12.4 


-27.2 


1.9 


fill 


N 


311 


7:21:00.52 


1:45:07.0 


17.09 


1.579 


-0.001 


11.1 


76.1 


2.1 


fill 


N 


389 


7:21:20.15 


1:45:57.2 


17.52 


1.679 




8.4 


3.2 


2.7 


fill 


N 


521 


7:20:52.79 


1:47:19.4 


16.96 


1.430 


+0.005 


11.7 


90.3 


2.0 


RGB 


Y 


538 


7:20:58.83 


1:47:31.2 


17.20 


1.577 


+0.003 


15.9 


97.4 


1.5 


RGB 


Y 


572 


7:20:54.89 


1:47:53.2 


16.66 


1.491 


-0.000 


10.5 


98.2 


2.2 


RGB 


Y a 


693 


7:20:57.23 


1:48:45.4 


17.14 


1.558 


+0.000 


5.2 


95.3 


4.1 


RGB 


ya 


706 


7:21:14.49 


1:48:54.0 


17.30 


1.461 


-0.047 


7.5 


2.4 


3.0 


fill 


N 
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FIG. 6. — Same as Figure 9 for Sal. 



FIG. 7. — Same as Figure 9 for Be39. 



been cataloged; Harris 1996), BH176 could provide impor- 
tant new constraints on Milky Way globular cluster formation 
scenarios. 

4.3. Be 29 

The RV determinations for Be 29 is within the errors (lcr) 
of other published findings (see table 11). For Be29, the mea- 
sured spectroscopic metallicity is Z=0.005 ([Fe/H] = -0.44; 
Carraro et al. 2004, hereafter CBVMPR). We explored 
Z=0.004 ([Fe/H] = -0.68), Z=0.008 ([Fe/H] = -0.38) and 



Z=0.019 ([Fe/H] = +0.0) Girardi et al. (2000) isochrones. 
Figure 10 shows the CMD for Be29 with the matched 
isochrones for the different adopted metallicities and for three 
different ages (3.2, 3.5 and 4.0) Gyrs. It is evident from the 
figure that both Z = 0.004 and Z = 0.019 isochrones do a poor 
job matching the slope of the RGB. The best matches are 
obtained for Z = 0.008, and we can see that the three ages 
seem to do similarly well matching the rather poorly defined 
MSTO, subgiant branch, and position of the RC. Ages outside 
this range fail to reproduce the MSTO and the position of the 
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FIG. 8. — Be29 and BH176 V/ photometry color-magnitude diagrams. 
Stars observed spectroscopically that were selected from the RGB in the core 
(< 2') are denoted with large open circles (o) for non-members and filled 
circles (•) are for RV members. Triangles (A) are stars observed outside 
the cluster radius or non-RGB selected stars, again filled triangles (A) denote 
stars selected as members. 



RC simultaneously. Together, the isochrone matches suggest 
a probable age range of 3.2-4.0 Gyrs, and E(V-I) ranging 
from 0.08 to 0.10. These values yield a heliocentric distance 
of 13.1-14.5 kpc. The resulting "revised" cluster parameters 
(Table 12) are roughly consistent with previous work. We find 
a similar age, reddening and distance for Be29 to that found 
by Tosi et al. (2004), though we find a lower reddening and 
thereby farther distance when compared with Kaluzny (1994, 
E(V-/)=0.10, d,Q = 10.5 kpc). Using the same isochrones for 
Be29 and Sal, we find that Be29 is the most distant known 
open cluster (Rcc = 21) in the Galaxy, as found by Tosi et 
al. (2004, see Table 12). The isochrone-fitted metallicity is 
in agreement, within our errors, to the spectroscopic, detailed 
chemical abundances of CBVMPR and Yong et al. (2005). 

However, Anthony-Twarog et al. (2005, hereafter AAT) 
pointed out in their comparison of NGC 2243 to Be29, that 
the findings of CBVMPR show that Be29 is somewhat a- 
enhanced. In Figure 11, we refit the CMD using the Padova 
a-enhanced isochrones (Salasnich et al. 2000). AAT matched 
a-enhanced isochrones to Be29 providing a reasonably good 
fit within previous determinations of the cluster redding and 



metallicity, although they do not have stars on the RGB higher 
than the RC. However as shown in Figure 11, our matches to 
a-enhanced isochrones are markedly different. While some 
of the differences can be attributed to the different isochrones 
used, the primary difference comes from trying to match our 
2 spectroscopically confirmed member stars above the RC. 
Matching the entire CMD with a-enhanced isochrones re- 
quires one to make Be29 have solar [Fe/H]. While this match 
is consistent with previous reddening determinations, it is 
very inconsistent with previously determined spectroscopic 
metallicities which range from [Fe/H] = -0.4 (CBVMPR) to 
[Fe/H] = -0.7 (Tosi et al. 2004). An additional method that 
comes near to matching the CMD, is to decrease the age to 
< 2.2 Gyr, however this match is inconsistent with the find- 
ing of AAT, who find that the age difference of Be29 and NGC 
2243 (Age = 3.8 ±0.2 Gyr) cannot be more than 0.5 Gyr. Ad- 
ditionally the reddening determined must be 50% higher that 
the highest previous estimate and double that found by AAT. 

We find that Be29 cannot be consistently fit with a- 
enhanced isochrones, therefore to account for the variable 
a element enhancements found in CBVMPR more complex 
models are needed. We refer the reader to §6.2 in Yong et 
al. (2005) for further discussion of the high resolution stud- 
ies. Using the same isochrones for Be29 and Sal (see be- 
low), we find that Be29 is the most distant known open cluster 
(R C c > 21) in the Galaxy, as found by Tosi et al. (2004, see 
Table 12). 

4.4. Sal 

Our Sal RV (V r , aU = 95.4 ±3.6 km s" 1 ) is smaller by ~ 
9 km s" 1 than the Sal velocity obtained by CBVMPR. We 
note that two of our Sal stars — 572 and 693 — were the 
only stars observed by CBVMPR for this cluster; CBVMPR 
find RVs of (104.4, 104.8) km s -1 for these two stars, respec- 
tively, while our survey finds (98.2, 95.3) km s -1 . Thus, while 
both studies find small RV differences between the two stars, 
the relatively large RV offset between studies persists and is 
apparently not related to a difference in the stellar samples. 
While CBVMPR does observe at much higher S/N than we 
do, the magnitude of the RV offset is significantly larger than 
the estimated uncertainty in the mean cluster motion for the 
two surveys. Moreover, we have found no net zero point off- 
set in checks of our standard star RVs by comparison to the 
IAU values to our measurements, as shown in Table 4. The 
same template stars were used for the RV determination of 
all of the clusters and no large offset is found in Be20, 29, 
39 or BH176. Thus, the origin of the systematic offset be- 
tween the surveys is not clear. However, since we observed 
the same stars as CBVMPR, we can determine that the four 
stars selected as members are indeed consistent with the work 
of CBVMPR. 

The sparse CMD of Sal, and lack of bright giants, lim- 
its our ability to achieve a reliable isochrone match, though 
our RV members provide crucial benchmarks at the magni- 
tude of the red clump. To improve the photometry depth and 
provide a better match, we used the photometry of Carraro 
& Baume (2003) within 2 ' of the cluster center and added 
our confirmed members. The spectroscopic metallicity for 
Sal gives a Z = 0.008 ([Fe/H] = -0.38; CBVMPR). We fit 
the ages 4, 4.5 and 5 Gyrs, with the fixed Z = 0.008 metallic- 
ity, are studied as shown in Figure 13. An age of 4.5-5 Gyr 
is derived, since the 4 Gyr isochrones do not match well the 
MSTO. This match results in anE(V-I) that ranges from 0.19 
to 0.24. The corresponding heliocentric distance for these val- 
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FIG. 9.— Girardi et al. (2000) isochrone matches for BH176. Metallicities of Z = 0.008 ([Fe/H] = -0.38), Z = 0.019 ([Fe/H] 
([Fe/H] = +0.2) are matched from top to bottom. Presumed cluster members, derived from spectroscopy are denoted with triangles (A). 



-0.0), and Z = 0.030 



ues ranges from 1 1.8 to 13.5 kpc (see Table 12). We find that 
we are in relative agreement about the cluster parameters with 
other photometric studies of Sal (Frinchaboy & Phelps 2002; 
Carraro & Baume 2003). Frinchaboy & Phelps (2002) find 
an older age (6.3 Gyr) based on the morphological age index 
(Janes & Phelps 1994), which is known to yield older age es- 
timates than those derived from isochrone matching. 

4.5. Cluster Membership 

With defined "RV-member" samples, we may examine the 
membership reliability of stars selected to be cluster mem- 
ber candidates by our photometric pre-screening. This in- 
formation is useful because a large number of our candi- 
dates have not been observed spectroscopically, but if we 
have confidence that the photometrically-selected candidates 
are likely to be members, they can potentially help con- 
strain the likely positions of the RGB and RC sequences in 
the cluster CMDs as well as point to the effectiveness of 
these kinds of selections for picking spectroscopic samples 
in future studies. The RV membership success rate among 
Washington+DD051 selected RGB stars with RVs is (100%, 
53%, 100%) for (Sal, Be39, and Be20) respectively. ERGB 



stars were added to our analysis for Be20 but none were 
found to be members. The membership success rate for 
the VI selected RGB stars was 71% for BH176 (60% for 
RGB+RC+ERGB) and 75% for Be29. In previous studies 
of dwarf spheroidal galaxies (e.g., Westfall et al. 2005; Ma- 
jewski et al. 2005), the Washington+DD<95 1 technique has 
been found to improve the selection reliability of star system 
RGB stars over CMD-only selection techniques by a signif- 
icant factor, even with samples selected over larger areas, to 
the low density extremities of the clusters. The present use 
of the Washington+DD<95 1 technique is not completely anal- 
ogous to these dSph studies, because here we have targeted 
any giant star candidates for spectroscopic follow-up without 
regard to CMD position (because in some cases the location 
of the RGB was not certain), whereas the dSph studies ap- 
plied both Washington+DD(95 1 selected of giant candidates 
and specific CMD location pre-filtering on the stars samples. 
Nevertheless, we do find that use of the Washington+DD(95 1 
filters has apparently improved our success rate relative to 
simply picking targets based on CMD position. For further 
discussion of the Washington+DD(95 1 selection technique the 
reader is directed to Majewski et al. (2000, 2005). 
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FIG. 10.— Girardi et al. (2000) isochrone matches for Be29. Metallicities of Z = 0.004 ([Fe/H] = -0.68), Z = 0.008 ([Fe/H] = -0.38), and Z = 0.019 
([Fe/H] = —0.0) are matched from top to bottom. Presumed cluster members, derived from spectroscopy are denoted with triangles (A). 



TABLE 1 1 
Bulk Cluster Radial Velocities 



Cluster 




BH176 


Berkeley 20 


Berkeley 29 


Berkeley 39 


Saurer 1 


# Members 


(RGB) 


4 


4 


6 


18 


4 


V r.RGB 


(km s~' ) 


9.1 


75.6 


27.1 


55.0 


95.4 


° V,, RGB 


(km s _1 ) 


4.6 


2.5 


3.8 


2.9 


3.6 


Vgsr,RGB 


(km s -1 ) 


-103.3 


-22.0 


-51.2 


-107.1 


-41.4 


# Members 


(all) 


9 


5 


11 


24 


4 


Vr,all 


(km s -1 ) 


11.2 


75.7 


28.4 


55.0 


95.4 


a V r ,aU 


(km s ) 


5.3 


2.4 


3.6 


2.5 


3.6 


Vgsy.a! 1 


(km s -1 ) 


-101.2 


-21.8 


-49.7 


-107.1 


-41.4 


Vr.pub 


(km s _1 ) 


85.0 a 


78.0 


24.7 


58.0 


104.0 


<?Vr,pub 


(km s _1 ) 


30 


5 


0.1 


2 


0.1 


Ret 




1 


2 


3 


2 


4 



References. — (1) F04; (2) Friel et al. (2002); (3) Yong et al. (2005); (4) CBVMPR. 
a RV based on one star, see §3.3 



5. DISCUSSION AND CONCLUSIONS disk as determined by old open clusters (e.g., Carraro et al. 

The clusters Be29 and Sal have potentially critical lever- l994 ^ Janes & Phel P s 1994 )' Even with lar § e statistical stud- 
age on radial, age and metallicity gradients in the Galactic most Galactlc radlal trends are hl g hl y dependent on the 

few open clusters in the extreme outer parts of the Milky Way 




FIG. 11.— Salasnich et al. (2000) Padova a-enhanced isochrone matches for Be29. Metallicities of Z = 0.008 ([Fe/H] = -0.38; bottom), and Z = 0.019 
([Fe/H] = —0.0; top) are matched. Presumed cluster members, derived from spectroscopy are denoted with triangles (A). 
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FIG. 12. — Girardi et al. (2000) Z = 0.008 ([Fe/H] = —0.38) isochrone matches for Sal, using the Carraro & Baume (2003) photometry. Presumed cluster 
members, derived from spectroscopy are denoted with triangles (A). 



(Carraro et al. 1994; Friel et al. 2002). As discussed in CB- 
VMPR and Yong et al. (2005), these two outer clusters deviate 
from the radial-metallicity trend found for the rest of the bulk 
open cluster population by Friel et al. (2002), which suggests 
either that the chemical properties of the outer disk maybe 
different from the rest of the disk, or that these clusters have 
nothing to do with the chemical evolution of the disk. We 
explore the latter scenario now. 

Friel (1995) has argued that the among the two most plau- 
sible models for the formation of the system of outer old open 
clusters, their creation as part of the normal secular evolution 
of the disk would require an unlikely "fine tuning" of forma- 
tion and destruction processes. On the other hand, formation 
of these clusters through accretion is a "natural mechanism" 



for forming these open clusters, which typically have large 
\Zgc\- We find that, when viewed together, various proper- 
ties of these old, outer open clusters lend circumstantial sup- 
port to their being part of a tidal debris system, and may well 
be "accretion products" from GASS. We compare the derived 
cluster galactocentric radial velocities (V gsr ) to the / vs. V gsr 
trend found for the GASS M giants and other clusters (Crane 
et al. 2003, F04) in Figure 14. Included are the analytical fits 
to the stream from Figure 2 of the analysis by Penarrubia et 
al. (2005) for comparison. The following equation is used to 
convert from the V r to V gsr trend, as used in F04 and Crane et 
al. (2003): 



V gsr = V r + 9 cos(b)cos(l) + 232 cos(b)sin(l) + 7 sin(b) (4) 
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TABLE 12 
Bulk Cluster Parameters 



Cluster 


Age(Gyr) 


d Q (kpc) 


Rgdkpc) 


E(V-I) 


[Fe/HL,, ec 


Vrikms' 1 ) 


V gsr (km s 1 ) 


GASS? 


BH176 


6.3 ± 1.0 


15.8 ± 0.5 


9.9 


0.70 




+11.2 


-101.2 


??? 


Berkeley 29 


3.7 ± 0.5 


13.4 ±0.4 


21.1 


0.09 


-0.44 


+28.4 


-49.7 


Yes 


Saurer 1 


4.5 ± 1.0 


13.1 ±0.4 


20.2 


0.22 


-0.38 


+95.4 


-41.4 


Yes 




FIG. 13. — The l-Vg Sr distribution of objects lying within 2.35 kpc of the GASS cluster plane and other clusters discussed in C03. Corrections to V gsr are 
described in §3.2. The solid lines are analytical orbits for the Monoceros/GASS stream taken from Figure 2 of Penarrubia et al. (2005). The hash marks represent 
a velocity dispersion of 30 km s -1 about the v g „ of the orbits, which approximately matches the velocity dispersion of GASS M giant stars from Crane et al. 
(2003). 
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BH176 




5 10 15 



Age (Gyr) 

FIG. 14. — AMR for all GASS clusters with spectroscopically determined 
metallicities, and RVs consistent with membership, with circles (•) denoting 
Be29 and Sal with values from this study. The solid line is the AMR for the 
Sgr dSph (Layden & Sarajedini 2000). The dotted line shows the metal poor 
limit of the basically flat AMR for the Milky Way old open cluster system 
(Friel 1995). The (x) denotes BH176 based on its photometric metallicity 
that does not fit the GASS clusters AMR. 



The RVs, when combined with the parameters in Table 11, 
yield a V gsr of (-50, -39, -101) km s" 1 for (Be29, Sal, and 
BH176) respectively. These RVs place both Sal and Be29 
within the GASS trend found in Crane et al. (2003) and, as 
shown in Figure 14, are generally consistent with the Mono- 
ceros models of Penarrubia et al. (2005). The new BH176 
RV is only, at best, marginally consistent with membership in 
GASS. As pointed out by Penarrubia et al. (2005) and Yong et 
al. (2005), additional information, especially proper motions, 
are needed to test whether these clusters are part of a coherent 
dynamical group consistent with the GASS/Monoceros sys- 
tem. 

Be 29 and Sal also follow the general age-metallicity trends 
suggested earlier for GASS (F04). We revisit (Figure 15) the 
GASS age-metallicity relation (AMR), now including only 
those clusters that are both spatially and dynamically con- 
sistent with GASS and for which spectroscopic metallicities 
have been derived: Be29, Sal, Tombaugh 2, Arp-Madore 2, 
Palomar 1, NGC2808 and NGC 5286. In addition BH176 is 
included (shown as x) to show that if the isochrone metal- 
licity is correct that it does not fit the AMR trend. Figure 
14 also includes the AMR trend derived for the Sgr system 
(Layden & Sarajedini 2000) and the metal-poor "limit" to 
the old open clusters of the Galaxy (Friel 1995) shown as 
the long dashed lines. The Sgr trend is typical of that ex- 
pected for an independently evolving, "closed-box" system 
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with protracted star formation, and it is clear that the pro- 
posed GASS cluster system more closely matches this Sgr 
AMR than the AMR of the nominal Milky Way open clus- 
ter system at smaller Galactocentric distances. We also com- 
pared to the metal poor envelope of the Carraro et al. (1998, 
small dashed line) AMR, which shows a "metal-poor spike" 
at ~ 3-4 Gyr. However, all of the three clusters that are more 
metal-poor than [Fe/H] = -0.6 have since had their metallic- 
ity re-evaluated and are found not to be as metal-poor. Still, it 
is interesting to consider that all of the clusters in this "spike" 
are all distant anticenter clusters, which requires one to ask 
whether this is due to a common formation, and if so what 
caused the "spike" 

We conclude that the RVs of the clusters Sal and Be29 
are consistent with GASS, which further points to an "ac- 
cretion" origin for these clusters, if one believes that GASS 
represents an accreting dwarf galaxy system. But whether 
this "accretion" origin is by the cluster being taken in by the 
Milky Way already intact or by it being formed in the process 
of accretion cannot be distinguish with current data. The 
detailed abundances of Yong et al. (2005) and CBVMPR lend 
weight to outer disk clusters being formed from an accretion 
event, since the abundance ratios (specifically [a/Fe]) are too 



high compared with current dwarf spheroidal galaxies. Be29 
and Sal may be associated with the thick disk which also 
has enhanced [a/Fe] (Venn et al. 2004) and could point to a 
common origin for both. 
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